Splashing onset in dense suspension droplets 
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We investigate the impact of droplets of dense suspensions onto a solid substrate. We show that 
a global hydrodynamic balance is unable to predict the splash onset and propose to replace it by 
an energy balance at the level of the particles in the suspension. We experimentally verify that the 
resulting, particle-based Weber number gives a reliable, particle size and density dependent splash 
onset criterion. We further show that the same argument also explains why in bimodal systems 
smaller particles are more likely to escape than larger ones. 



Splashing of liquid droplets upon impact on a solid sur- 
face has been investigated for over a century [1 . A long- 
lasting pursuit of finding the parameters governing the 
splashing transition identified factors ranging from the 
droplet impact speed to the roughness of the impacted 
surface to the pressure of the surrounding gas pHT2]. 

More recently, there has also been a growing interest in 
what happens to the spreading and splashing if particles 
are added to the liquid or when a non-newtonian fluid is 
used p!3HT8] . On micron scales, Zr02 suspensions have 
been used in studies aiming to optimize ink-jet print- 
ing applications [HHSS] , and on truly macroscopic scales 
there has been the development of 3D printers that dis- 
pense cement slurry [24H26] . In all of these situations, 
an important concern is to prevent splashing, and par- 
ticles from escaping, when droplets hit a surface. How- 
ever, the question of when and why particles are ejected 
has remained unsettled, and existing experimental stud- 
ies mostly focus on dilute suspensions. 

Current models for suspension drop impact associate 
the onset of splashing with the condition that K = 
We]^'^ Re]/^ exceeds a critical value i^o, which has been 
the traditional criterion for pure liquid splashing on dry 
surfaces at atmospheric pressure in a regime independent 
of surface roughness [4^ 27^ 28 . Here the droplet-based 
Weber and Reynolds numbers are defined as We^ = 
pi^dU'^/o' and Red = piVdU/fi, with rd the droplet radius, 
U the impact velocity, and pi, a and p the liquid density, 
surface tension and dynamic viscosity respectively. 

In these models, the addition of particles to a pure 
liquid has been captured by replacing p with an effec- 
tive viscosity /ig that increases with the packing frac- 
tion of solid particles [14| [29H33] . This predicts that K 
should decrease as particles are added: a droplet of a 
pure liquid that would splash under certain conditions 
should not splash after adding enough particles. To our 
knowledge, there exists no systematic study that confirms 
this prediction. Surprisingly, Nicolas [14 observed that 
adding particles, instead, lowered the splashing thresh- 
old Kq. This contradicts the energy balance based on 
global, droplet-scale properties and clearly suggests that 
the splashing onset is affected by particle- level aspects. 

To investigate how impact energy is transferred to the 



particles that are ejected, we depart from the dilute limit 
described above, but instead start at the other side of the 
spectrum, i.e., we will focus here on the limit of dense 
granular suspensions with volume packing fractions (j) = 
0.62 ±0.03, where the discrepancy with the droplet-scale 
splash onset criterion is most pronounced. 

There is, however, a clear difference between the splash 
of a pure liquid and a suspension. In the liquid case, 
the size of the eject a depends on either the destabiliza- 
tion of a thin liquid sheet I34H39] or, in the case of 
prompt splashing, an instability at a moving contact 
line [5, 11, 40 . At splash onset in a suspension, on the 
other hand, the ejecta are the solid particles (see Fig.jl]), 
which implies a built-in length scale. This length scale 
was not taken into account in the energy balance lead- 
ing to K, and not considered by Refs. JHH |23]. We 
will now evaluate the energy balance at the particle level, 
which we will then compare to our experimental obser- 
vations. 

Energy balance -Suifdice tension keeps particles inside 
because an escaping particle involves an increase of the 
surface area. A particle can thus escape if it has enough 
kinetic energy to overcome surface tension. In compar- 
ing the surface energy associated with escaping particles 
to the kinetic energy of particles we assume that during 
the impact of a suspension drop the velocities are high 
enough to neglect dewetting of particles [4T1 . This means 
that an escaped particle still has a (thin) layer of liquid 
around it and the increase in surface area will scale with 
the surface of a particle: 

Esurf ^ 47rr^cr, (1) 

where is the particle radius and a is the surface tension 
of the liquid. The kinetic energy of a particle is given by 
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Ekin = -TTppr^ul, (2) 

where pp is the specific density of the particles and Up 
the velocity of a particle. The velocity Up is a result 
of collisions between neighboring particles (see Fig. [2|, 
which converts the vertical velocity into a horizontal one. 
Based on momentum conservation we expect that the 
velocity Up of a particle sitting on the outer surface of a 
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FIG. 1. Still images right after impact of Zr02 particles in water onto a smooth glass plate, shown for different droplet radii 
rd (see text above the images) and particle radii: (a-i) rp = 362 /xm, (j-l) Vp — 138 /xm, (m-s) rp = 78 /xm. The time between 
impact and ejection of the first particles ranges from 0.6 to 3.5 ms. Images (a-o) are organized in vertical columns, with drop 
impact speed and therefore K increasing from top to bottom in order to bracket the onset of splashing, defined here as the 
ejection of individual particles. For the K values listed in the top/middle/bottom rows we never/sometimes/always observed 
particle ejection. Here K — We^J.'^ Re^^^ ^ where we used the effective viscosity \ie given by Krieger and Dougherty [31] in 
the Reynolds number. The blurred background in (a) and (6) is the out-of-focus image of the syringe. We note that (o) has 
been thresholded and dilated in order visualize the ejected particles that would otherwise be invisible due to their small size, 
(p,^) side and bottom view of a droplet that does not splash (Wep — 12), (r,s) side and bottom view of a splashing droplet 
[Wtp — 26). The scale bar in the images is 5 mm, images (a-6) all have the same scale. 




FIG. 2. Sketch of the ejection mechanism upon impact. 



droplet will be similar to the impact velocity of the drop 
U . The ratio of the kinetic and surface energy then is 
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Here and throughout the rest of this Letter, Wcp is a 
particle-based Weber number defined by the density and 
the radius of the particles, the impact velocity of the 
suspension droplet and the surface tension of the liquid. 
Fig.[TJp-s) shows an example at We^ = 12 where the sur- 
face tension is strong enough to keep all particles inside 
the suspension and at Wcp = 26 we see that particles are 
ejected in all directions. 

Experiments -We prepared suspensions of demineral- 
ized water with Zr02 (Glenn Mills) and soda-lime glass 
(Mo-Sci) particles inside a syringe. All particles are 
spherical, with standard deviations from their mean size 
of 5 to 8% for the Zr02 and 15 to 20% for the glass 



beads. The packing fractions were determined by mea- 
suring the mass of a suspension drop, letting the water 
evaporate on a hot plate, and then measuring the mass 
of the dry particles. All packing fractions were between 
59 and 65%. The packing fractions are probably slightly 
overestimated due to the finite sample size, because the 
amount of liquid depends on the shape of the menisci be- 
tween the particles sitting at the surface of the droplet. 
In the current study we focus on the changes to inviscid 
liquid splashing introduced by the particles and do not 
explore the role of additional viscous dissipation from the 
suspending liquid [42] . 

The drops were formed by very slowly pushing the sus- 
pension out of a syringe using a syringe pump. We cut 
off the narrow end of the syringe, resulting in a sim- 
ple, straight cylindrical nozzle. As gravity pulls the sus- 
pension down, a pinch-off will occur g3[ |44], resulting 
in highly reproducible suspension drops. We varied the 
impact speed by adjusting the release height, and drop 
sizes by changing the nozzle size. During the extrusion 
of the suspension there is minimal deformation of the 
droplet, which results in a droplet radius equal to the 
nozzle radius . The experiments were recorded with a 
Phantom V12 high speed camera at frame rates of 6,200- 
10,000 fps with a 105 mm Micro-Nikkor lens, resulting 
in a typical resolution of 20 — 50 /im/pixel. Side views 
were used to determine the impact speed, and a bottom 
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view was used to observe and track particles that are 
ejected after impact. Typical bottom view images are 
shown in Fig.[l](a-o). We never observed the ejection of 
liquid droplets in our experiments. 

To determine the transition velocity /7* above which 
particles are ejected from the droplet, we determine by 
repeated experiments the lower bound of the impact ve- 
locity at which we never see a splash and an upper bound 
at which we always see ejected particles. Because we are 
able to distinguish individual particles, we define a non- 
splashing experiment when not a single particle leaves 
the suspension. A particle has left the suspension when 
there is no liquid bridge connecting it to the other par- 
ticles. The lower and upper bounds are represented by 
the error bars in Fig. [3j The experimental determination 
of [/* for the case of Vp = 78 /im and Vp = 362 jam Zr02 
particles is shown in Fig. [4ja), where Ns is number of 
times we observe a splash and N is the number of times 
we repeat one impact speed U (typically 10 times). 

Results on unimodal suspensions-Fig. [l] verifies that 
a global criterion, K > Kq, does not capture the ob- 
served behavior. Within these examples there is about 
a factor 5 difference in i^o-values for droplets comprised 
of the same liquid and similar packing fraction of Zr02 
beads. Note that in Fig. [ija-^) the particles also all have 
the same size. Fig. [3] shows the infiuence of the particle 
size on the splash transition. In all cases, the transition 
to splashing happens at the same value of Wep ~ 14. 
This is consistent with Eq. (|3| , where we argued that the 
particle-based Weber number is the relevant parameter 
for the splash onset [46]. Possible non-Newtonian effects 
or an effective viscosity seem to play no role here. Com- 
paring the black, cyan, and magenta lines in Fig. ^a) 
-corresponding to the cluster of data points in the far 
right of Fig. [5]- we see that a factor of over four in the 
droplet radius has no significant infiuence on the splash 
onset - in strong contrast with using the global droplet 
size as relevant length scale. 

The measured velocity Up of an ejected particle is al- 
ways similar to the impact speed, confirming our estimate 
used in Eq. (|3|. The inset in Fig. [3] shows a typical veloc- 
ity distribution for 78 /im particles at Wep = 20, where 
the mean velocity Up of the particles after they are ejected 
from the suspension is slightly higher than the impact ve- 
locity U. For the larger Zr02 particles that we tested Up 
was smaller than U and the ratio Up/U decreases with 
particle size, but all measured mean velocities are in the 
range 0.5 < Up/U < 1.2. 

Results on bimodal suspensions-Fiom Fig. ^a) it is 
clear that a suspension with particles of radius 362 jam 
always splashes above an impact speed of about 1 m/s. 
On the other hand, a suspension with particles of radius 
78 /im never splashes below an impact speed of about 
1.3 m/s. So what happens if we make a bimodal suspen- 
sion of these two particles types and impact at a speed 
between these two limits? We find that the splashing 




FIG. 3. Splash onset velocity U* as a function of the product 
of particle radius rp and particle density pp. The cluster of 7 
data points at the far right correspond to 3 different drop sizes 
(see also Figs.[l]|4| and substrates with 4 different roughnesses 
(see main text), all other data points are for = 2.3 mm on 
a smooth substrate. The solid curve gives the onset Weber 
number Wcp* calculated from a best fit to all experimental 
data. The dashed lines represent the upper and lower bounds 
for Wcp* , corresponding to one standard deviation. The error 
bars of the individual onsets are obtained by determining the 
upper and lower bound for [/* corresponding to the spread in 
the experimental data. Inset: Histogram of the velocity Up of 
ejected particles at Wcp = 20, for rp = 78 /xm. 



behavior is inverted: The larger particles remain inside, 
while the smaller particles get ejected. In Fig. ^b) we 
determine the splashing behavior for the bimodal suspen- 
sion the same way as we did for unimodal suspensions 
in Fig. ^a). Clearly, the transition curves switch their 
position when going from the unimodal to the bimodal 
suspension. 

At first sight this result seems counterintuitive: surface 
tension should be more effective in keeping small particles 
inside the droplet than large particles. To explain this we 
take a closer look at how the particles obtain their ve- 
locity upon impact. Just before the moment of impact 
all particles have the same velocity U directed down- 
wards. Then, upon impact, the vertical velocity needs to 
be transferred to horizontal velocity. This transfer is a 
result of collisions between particles. The velocities that 
are obtained from these collisions obey momentum con- 
servation. Collisions between particles of the same size 
will result in similar velocities - which explains why Up 
scales with U. Collisions between large and small parti- 
cles can, however, result in much larger velocities for the 
small particles [471 |48j , which explains why the transi- 
tion is at a lower impact speed. Conversely, small par- 
ticles can only give little velocity to large particles. The 
presence of small particles also reduces the chance for 
direct collisions between large particles, which explains 
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FIG. 4. Comparison between unimodal and bimodal Zr02 
suspensions. All data for droplets made with nozzle radius 
Td — 2.3 mm, except for the data shown in magenta and 
cyan in (a) for which Vd — 1.0 mm and Vd — 4.3 mm, re- 
spectively, (a) Splash onset for unimodal suspensions of large 
(rp = 362 /im, black, cyan, and magenta open symbols) and 
small {vp — 78 /xm, red dots) particles, {h) Splash onset for 
bimodal suspensions with volume ratio of small to large par- 
ticles about 1:1. Two distinct onsets exist; one where only 
small particles escape and one where both small and large 
particles escape. Data symbols and colors as in (a). All insets: 
Examples of the unimodal and bimodal suspension droplets 
corresponding to the data presented in this figure, imaged 
from below just before impact, (c) Histogram of the relative 
momentum change for large particles in unimodal (cyan) and 
bimodal (black) suspensions at [/ = 2 m/s. (d) Same as (c), 
but for the small particles. 



the increased transition impact speed for large particles 
in bimodal suspensions. 

In order to extract the change in momentum due to 
the impact, we perform experiments for all three suspen- 
sion configurations at an impact speed of 2 m/s where 
we always find ejected particles: (i) unimodal suspen- 
sions with Vp = 362 /im, (ii) unimodal suspensions with 
Vrp = 78 /im, and (iii) bimodal suspensions consisting 
of the same two different particles sizes. We calculate 
the gain in momentum due to the collisions as follows: 
Because all these experiments are above the splashing 
transition we can directly determine the velocities of the 
ejected particles and thus their kinetic energy. We know 
that during the ejection process the particle has trans- 
ferred part of its kinetic energy to the surface energy 
given by Eq. ([T]). Adding this surface energy to the mea- 
sured kinetic energy gives us the kinetic energy Ein of 
the particle just before it was ejected, which corresponds 
to the momentum = (2m£^^n)^/^, with m the mass 
of an individual particle. Comparing this to the vertical 
momentum of the particle before the moment of impact. 



Po = mU, gives us the relative change in momentum 
Pin/Po due to the collisions between particles. 

In Fig. [4]^c-(i) we show the ratio Pin/po for the three 
experiments mentioned above. Every impact speed was 
repeated 10 times and we count the number of ejected 
particles rip that are within a specific range of momenta. 
Fig. Qc) clearly shows that for large particles going 
from unimodal suspensions (cyan) to bimodal suspen- 
sions (black), the probability of finding particles with a 
momentum in the order of po decreases by at least an or- 
der of magnitude. Note that for the bimodal suspensions 
we found for all large particles Pin/po ^ 0.75. For small 
particles on the other hand (Fig. |4j(i)), we observe that 
the probability of finding small particles with a momen- 
tum ratio > 1 is increased by more than a factor 5 when 
changing from a unimodal to a bimodal suspension. 

Since collisions between particles are responsible for 
driving the onset of splashing, and less so the interaction 
of particles with the substrate, we expect that surface 
roughness will play a much smaller role than for pure 
liquids. To check this we have performed experiments 
on rough substrates with roughness length scale i rang- 
ing from £> Tp to i <^ Tp. The results are included in 
Fig. [3j where it is clear that in all cases there is no differ- 
ence compared to the impact on a smooth surface. The 
same independence of the splashing onset holds for exper- 
iments performed at a reduced pressure of P ~ 10 kPa. 

Conclusions and outlook -These experimental results 
demonstrate that the relevant parameter to quantify the 
splash onset is a Weber number calculated at the particle 
level. This is in contrast to earlier proposals for splash 
onset criteria in suspensions. Local interactions between 
particles at the edge of the droplet are responsible for the 
ejection of particles upon impact. This explains why the 
effective viscosity, which acts on a global droplet level, 
does not prevent the suspension from splashing. Our 
observations give rise to the question at which packing 
fraction the global droplet description breaks down and 
when the dense limit, investigated here, takes over. The 
same local interactions are responsible for an inversion 
of the splash onset for bimodal suspensions. Since the 
momentum transfer in dense suspensions is dominated 
by collisions between particles, it might be possible to 
further tune the splash onset via particle characteristics 
that affect collisions, including particle shape, restitu- 
tion coefficient, or friction. However, our findings already 
demonstrate that splash onset in dense suspensions be- 
haves qualitatively different from predictions based on 
pure liquids. The typically used K parameter does not 
delineate the splash onset properly and neither substrate 
roughness nor ambient air pressure play the same sig- 
nificant role. Instead, a particle-based, critical Weber 
number Wep* describes the onset well. 
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